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Abstract 

We study the high-energy limit of tree-level string production amplitudes 
from decaying D-branes in bosonic string theory, interpreting the vertex oper- 
ators as external charges interacting with a Coulomb gas corresponding to the 
rolling tachyon background, and performing an electrostatic analysis. In partic- 
ular, we consider two open string - one closed string amplitudes and four open 
string amplitudes, and calculate explicit formulas for the leading exponential 
behavior. 
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1 Introduction 



A famous feature of string theory is the soft behavior of scattering amplitudes in the 
high-energy limit. Remarkably, at each order (worldsheet genus) in the perturbation 
theory, the dominant saddle point contribution has a simple electrostatic interpre- 
tation - the exponent can be identified as the electrostatic energy of point charges 
at the equilibrium jT]^3] . Recently, the electrostatic approach was applied to string 
scattering from decaying D-branes |4j[5]. In the original work jT]^3] , there was no 
background charge in addition to the point charges, due to momentum conservation. 
A new ingredient for strings scattering from a decaying D-brane is the condensing 
tachyon on the brane, which provides a background charge distribution that interacts 
with the point charges. More precisely, the background consists of a Coulomb gas of 
unit charges at finite temperature and (imaginary) chemical potential. This electro- 
static interpretation is valid at all energies. Exact calculations are very complicated, 
but simplifications can be found in the high-energy limit. In a previous work, we 
applied this approach to study closed string pair production from a decaying D-brane 
at disk amplitude level [5]. 

In this paper, we study further examples of the electrostatic approach to decay 
amplitude calculations. We start by a detailed general description of the approach. 
For a uniform discussion, we review previous calculations before applying the method 
to new cases. We verify that the approach predicts the known high energy behavior 
of the bulk-boundary amplitude [6j, and show that the amplitude does not diverge 
at high energies in the kinematically allowed region. We also present an improved 
analysis of the ( "antipodal" ) two-point function, related to open string pair produc- 
tion by the decaying brane. Then we discuss n-point functions with n > 2, and the 
associated string production amplitudes^ As a completely new result, we apply our 
method to calculate the high-energy limit of a three-point amplitude involving a pair 
of open strings and a closed string. For stable D-branes, such amplitudes are, e.g., 
relevant for emission from or absortion by a charged black hole [?] . Finally, we extend 
previous work on n-point boundary functions [8] and study the high-energy asymp- 
totic behavior of the amplitude for production of four open strings, before ending 
with conclusions and outlook. Technical details are presented in two appendices. 



2 Amplitudes from electrostatics approach 
2.1 Preliminaries 

We review quickly the structure of the disk amplitudes for string scattering from a 
decaying D-brane in bosonic string theory, with the "half S-brane" rolling tachyon 
background [9,10 . The correlation functions are path integrals 



t- n j2 I n 



w a ,w a )e 5Sbdry 



a=l \a=l I f rcc 

1 In this paper we consider only tachyon vertex operators. 
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where the rolling tachyon boundary deformation representing the half S-brane decay 
mode for the D-brane is 

SS hdly = X f dte x °^. (2) 



The calculational complications arising from rolling tachyon background are similar 
for all open and closed string vertex operators, when they are represented in timelike 
gauge 11 : 

V(k) = V sp (k)e- lujX0+it * , (3) 

— * — * 

where V sp (k) involves only space directions X. The contributions from the con- 
tractions from the spacelike sector are similar to those in scattering from stable (or 
non-deformed unstable) D-branes, as they do not involve the timelike background ([2]). 
Our focus is in developing techniques for calculating the contribution from the non- 
trivial timelike sector. Therefore, for our purposes it is sufficient to focus on scattering 
amplitudes which involve only closed and open string tachyon vertex operators 

V{w a ,w a ) = e ik °- X( - W ^ . (4) 

Closed string vertices are placed in the interior of the unit disk \w a \ < 1, whereas the 
open string vertices lie at the boundary w a = e tTa , r a = . . . 2tt. We will consider 
open-closed n-point amplitudes, with n = n c +n where n c (n ) is the number of closed 
(open) strings. We adopt a notation £ a = — iu a and break up the spatial momentum 
to parallel and perpendicular directions to the unstable Dp-brane: k a = (kt, k^). On- 
shell conditions for the bosonic closed string tachyons are k 2 = (,l + (k a ) 2 = —4(N a —l) 
and for the open string tachyons they read k 2 a = ^ + (kl) 2 = —(N a — 1) |6|. The 
overall spatial parallel momentum will be conserved: $2o=i ki = 0. 

The worldsheet correlation functions can be evaluated by first isolating the zero 
modes from the oscillators, = x M + X ,fl , and then expanding the boundary defor- 
mation into a power series in A. This yields 

N n c , 2 „ n 



V AT n * J 1 « nr, I 



dr a 

. . , - 2tt 

N=0 a=l a=n c +l 



N ,, IN 

?X'°{ti) TT Jk a -X'{w a ,w a ) 



dti I T—r \ ■'■),. 




(5) 



where we introduced z = 2n\e x and fixed the indexing of the vertex operators such 
that closed strings have smaller values of a. 
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After a straightforward calculation (see Appendix [AJ , we can write ^ as 

A nc+no = J dx°d p xh l ^ k ^A nc+no (x°) (6) 

oo 

A nc+no (x°) = J2(- z ) Nl ^ + n ({w a ,k a };N) (7) 



N=0 

oo „ n c , 2 



£(->Vn^ n if n 

Af=0 a=l a=n c +l l<a<b<n 

n, 

i r i _/../. _l /'— . /'— 

l<a<fe<n c a, 6=1 

xZ„ c+ „, o ({w a , /c a }; AT) , 



Fa - 



|"| \w a -W b \ k *- k l +k "-K Yl |1 - W a W b \^ k "' k b k "' k b"> 



where 



Z nc+no ({w a ,k a };N) = ^ J] K-^| 2 ^ 6 (9) 

l<a<fe<n 

r N m n 

i=l l<i<j<N a=l 

2.2 Applying the electrostatic approach 

In general, ^ is too complicated to calculate, so we look for an approximation 
scheme in a special limit. As suggested by the notation, Z nc+rio can be interpreted as 
a Coulomb gas partition functior]^]at the inverse temperature (3 = 2. We can evaluate 
the partition function in a saddle point approximation corresponding to electrostatic 
equilibrium, at large N and £ a (with the ratios £ a /iV fixed). The method is discussed 
in detail in a companion paper [4] . The leading term of log Z„ c+no can be found by 
evaluating the electrostatic energy of the Coulomb gas, unless n c + n Q is too large. 
The next-to- leading term was also analyzed in [2] . In this article we shall continue the 
analysis of [i] by discussing in detail how the high energy limit of string amplitudes 
arises from the electrostatic approach. String production at large £ a = — iu a is a 
physically interesting regime, since we expect the unstable D-brane to mainly decay 
to very massive string modes 
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First we note that the correlator ~ Z nc+Tlo needs to be summed over N and 
integrated over x° to obtain the string scattering amplitude, 

/oo 
J2(-z) N Z nc+no ({ZaY,N) , (10) 
JV=0 

where z = 2n\e x ° and we dropped the complicated proportionality factor for clarity. 
As discussed in [8], under certain conditions the sum and the integral can be carried 



2 We have explored many other aspects of this connection for decaying D-branes in bosonic and 



superstring theories in 12 ■ 16 
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out exactly, 



N=0 



(11) 

This result was originally found in [17] in the context of Liouville field theory. Z nc+no 
is in general unknown, but simplifies in the large N limit (N,£ a — > oo with N/£ a 
fixed) g]. Since A" = — ^ a £ a in , the large N limit corresponds to the high 



energy (u a = i£ a ) limit for the string production amplitudes. In particular, the 
leading term gives the saddle-point approximation 

7r(27rA) _ ^^ a f _ r i 

A nc+ n = . f v n (dwdwdT)~f(w,w,T)exp -2£\ N= _ T * , (12) 

where £ is the electrostatic energy of the corresponding Coulomb gas configuration 
in the continuum limit, and can be computed explicitly |4|. Above j (dwdwdr) de- 
notes the w a and r a integrals, and j(w, w, r) denotes the contribution from spacelike 
contractions in Their explicit form depends on the particular amplitude. Note 
that in general also £ is w a , r a dependent. If possible, we perform the integrals explic- 
itly, but in general we follow the approach offl] and replace w a by their electrostatic 
equilibrium values (for details, see Appendix |B| and the examples in Sections [3] and [4]). 

The use of (11) requires an analytic continuation of Z ric+rio (N) in the parameter 
A" to noninteger N, which is subject to some constraints. The partition function 
should be analytic for Re A" > 0, and it should not grow exponentially as A" — > oo, 
when |argA^| < 7r/2. The asymptotic behavior of Z nc+rio (with £ a fixed) may be 
analyzed using random matrix theory techniques. Using a random matrix theory 
interpretation, Z nc+no describes an expectation value of a periodic function in the 
circular ensemble of U(N) random matrices, CUE(A^). The expectation value can 
in turn be converted to a Toeplitz determinant of the Fourier coefficients of periodic 
funtion (see (6||18||19) for more discussion), the advantage then is that the asymptotic 



behavior of the determinant at iV oo simplifies, and is described by the Szego 
formula and its generalization. Using the result of (20l , we find the power-law behavior 



Z nc+no {N)~ (13) 



7V->oo 



(see [5] for a discussion with n c = 0). This result holds also for the leading term in the 
A^ — > 00 limit when N/£ a is fixed, an explicit discussion for the boundary one-point 
partition function is given in j2]. 

The above analysis needs to be modified if the large A^ limit of the partition 
function is not analytic in N. As explained in |4|[5], nonanalytic behavior may be 
associated with suddenly appearing or disappearing gaps^] in the continuous charge 
distribution of the Coulomb gas picture, created in the vicinity of the external charges 
at some large value of N. This may happen in the presence of bulk charges: a charge 
near the boundary of the disk generates a gap, which disappears at sufficiently large 

3 A gap may also exist first and then disappear, or many gaps may be created or join together. 



5 



N [5). Similar phenomenon may occur in higher order boundary amplitudes: two 
nearby charges will create only one gap for small N, which separates into two gaps 
as N increases. In these cases, as suggested in [5], one can first integrate over the 
locations of the charges. The integrated partition function is analytic in N, and one 
can proceed with the analytic continuation]^] 

In order to calculate the electrostatic energy, we need to first assume that £ a are 
real valued, and then continue analytically to physical energies u a = in the end. 
We justify this by noting that Z nc+no is analytic for Re£ a > —1/2, as seen from Q j^] 

In summary, we have developed a prescription for calculating the high energy 
approximation of the string scattering amplitudes, containing the following steps: 

1. Find the series expansion for the scattering amplitude (see Eqs. (J6|)- (J9]) ) and 
identify the Coulomb gas partition function. 

2. Solve the electrostatic potential problem to find S, and the leading partition 
function log Z ~ — 28 at (3 = 2. 

3. If necessary, integrate the leading result over the positions of the external 
charges (the unfixed modular parameters of the string vertex operators). 

4. Do the summation over N and integration over time by continuing analytically 
iV to the total external charge N — > — £ a in the result for log Z, as shown 



in Eq. (11). 



5. Continue analytically to physical energies u a = i£ a . 

Finally, let us comment on the precision of the obtained approximation. We only 
used the leading term of logZ nc+rio in the large N limit, which was shown to be 
0(N 2 ) in pi. Since we fixed N = -J2 a ^, our result Q is the leading O (w 2 ) 



term of log A nc+no in the limit of large u = ^2 oo a with the ratios oo a / 00b fixed. We 
also calculated the next-to-leading O (N) term of \ogZ nc+rio in [2]. However, this 
term vanishes at = 2. Therefore, the first nontrivial corrections to our high-energy 
approximation are obtained from the O (N°) term of log Z nc+Uo , resulting in O (u>°) 
(possibly logarithmic) corrections to \ogA nc+rio . 

2.3 Electrostatic energies 

In the remainder of the paper we study the easiest string scattering amplitudes follow- 
ing the above systematic prescription. In the second step, we need the electrostatic 



potential energy £, as in Eq. (12). We have already studied this problem in 4 5 



4 Another method which typically avoids non-analyticities is to fix the positions by the equilibrium 
equations, which follow from minimizing the electrostatic energy and will be discussed in Appendix 
B. 

5 There may be one caveat: the amplitude contains an integration over the unfixed moduli pa- 
rameters, which might in principle cause problems in the case of higher-point amplitudes, since the 
integrals are not necessarily well defined for all Re£ a > 0. This needs to be checked case by case. 
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and summarize the results here. First, the energy with one bulk charge £ = £1 at 
w = w\ = r (0 < r < 1) reads [HI 



gCbulk) = e ( rc _ r ^ ]og ( 1 _ r 2} ( M ) 



+9(r — r c ) 



(iV + 2£) 2 l + X N 2 l-x 



4 ° 1 + 5(r) 4 ° 1 - 5(r) 



+ ^rlo. 



2 & (l + 5(r)) 



where 



N 1 - r £ 

' *(0 = ^— ; x = tt~ ^ • (15) 



iV + 2£ ' v/ 1 + r ' * N + £ 
Several configurations with charges on the boundary of the disk were solved in In 
the simplest case there is one boundary charge £ on the boundary, giving 

£ (lpt) = ^[-2F(£) + F(2£)-F(iV)-F(iV + 2£) + 2F(iV + £)] , (16) 

where ^ ^ 

F(x) = -xf(x) = -x 2 \ogx . (17) 

We will also use the result for the configuration where two boundary charges £i and 
£2 lie at exactly antipodal points on the unit circle. Then 



£ (2pt) = {i ? (iV + £ 1 + £ 2 )-F(£ 1 + £ 2 ) 



_i [F(iV + 2£ x + 2£ 2 ) + F(N + 2£ x ) + F(N + 2£ 2 ) + F(N)] 
+ - A [F(26 + 2£ 2 ) + F(2£0 + F(2£ 2 )] - ^2 log 2 

As an example of a more complicated situation we consider a symmetric four-point 
case where two particles with equal charges £1 are, say, at r = axgw = and at 
r = 7r, and two additional particles with charges £ 2 are at r = 7r/2 and at r = 3tt/2. 
This configuration gives 

£(4 P t) = 2F(iV/2 + £x + £ 2 ) - 2F(£x + £ 2 ) 

- l - [F(N/2 + 2£x + 2£ 2 ) + F(N/2 + 2£ x ) + F(N/2 + 2£ 2 ) + F{N/2)\ 

+ l - [F(2£ x + 2£ 2 ) +F(2£x) + F(2£ 2 )] - +£ 2 ) 2 log2 . (19) 

3 Two-point amplitudes 
3.1 Bulk-bulk amplitude 

The high-energy production amplitude of two closed strings from a decaying D-brane 
was derived in [5] by using the framework of Section [2] with the result of Eq. (14). 
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For completeness, we review the results here. We fix the charge £ 2 at w 2 = 0, and 
the charge £1 at < W\ = r < 1 by using conformal symmetry. In this case Eq. (12) 
becomes 



4)+ 2 (a,&) 



7r(27rA) 



-6-6 



sin7r(£i + £ 2 
where we restored the proportionality factor 



drk{r) exp [-2£ (bulk) ] 



k(r) 



fci-fe 2 , 



1 — r 



2^(( fe i') 2 -(^) : 



(20) 



(21) 



For the special case where the strings have equal energies, U\ = uj 2 = w, the result 
simplifies to 



where 



ine 1 ™ (27rA) 
2 sinh 2nu 



2iui 



^ rQ- 1 ) , r(f + i) 
r(t±«_2) r(^) 



J [ 2~ l 



(22) 



s 
t 

u 



ki ■ k 2 

2(M') 2 = 2(4') 2 

(ki - kif 



_L\2 



(kty - (kt) 



(23) 



In the definitions of the parameters we used the on-shell conditions k\ = —u 2 + k 2 = 
— 4(iV a — 1) for the tachyonic states N a = 0. Notice that the results for the leading 
high energy asymptotics (see [5] for an extensive analysis) remain valid also for higher 
iV a 's as long as they are much smaller than the (squared) energy scale. 



3.2 Bulk-boundary amplitude 

Let us then discuss the asymptotics of the bulk-boundary scattering amplitude, and 
verify that our result mathces the exact one §. We follow §[Il)[l8] and use con- 
formal symmetry to place the bulk operator (charge £ c ) to the origin, rather than 
integrating over its position as suggested in (JsJ) . Then the bulk charge decouples from 
the Coulomb gas calculation, and we may use the results with one boundary charge 
(£ D ) for the partition function. According to Eq. (12), analytic continuation of (16) 
£o at 



to N 



Al+i(6o Co 



2 gives 

7r(27rA)-k 



sin vr (f c + £ 

Setting here £ a = — iu a we find the asymptotics 

in (27rA) i( " c+Wo) 



exp[-2 £ (lpt) ] 



t=to,N=-t c -to 



(24) 



A 1+1 (u ,uj c 



sinh 7r (oj c + lo 
1 
2 



exp 



u 2 Q log ■ 



u 2 c log 



u ) 2 loj. 



u 2 (21og2±ivr) + C(logw) 



(25) 
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where we restored the expected size of the next-to-leading order correction. This 
indeed matches with the asymptotics of the exact amplitude [6] up to the branch 



choice of the logarithm (± in the phase factor, on the last line in (25)) which is hard 



to obtain from the electrostatic approach. Notice, however, that the absolute value 
of the amplitude is independent of the branch. 

Let us make one comment about this result. After using momentum conservation 
parallel to the D-brane, 26-momenta of the strings become 



k c = {u c ,h\k^) 
k Q = (u , -k\0) . 



(26) 



At high energy, and for low-lying excitations (N a <C co 2 ), the mass-shell conditions 
-cu 2 + k 2 c = -4(iV c - 1) and -u\ + k 2 = -(N Q - 1) give 



UJ r 



UJ 



k\ 



(27) 



so asymptotically u c > u Q . The leading term in (25) can be written as 

u 



u a log 



U c + L0 o 



+ u c log 



00, 



s2, W, 

oj ) log 



L0 C + L0 o 



^"+2u 2 log 2 



(2* 



i 



a 2 log a + (l-a) 2 log(l-a) - - (l-2ay log(l-2a) + 2a 2 log 2 



where a = uj / (u c +uj ). In the kinematically allowed region < a < 1/2 the function 



in the square brackets in (28) is negative, and it vanishes at the endpoints a = 0, 1/2. 



Thus the amplitude vanishes for high energies in the kinematically allowed region 
as A ~ e ~<uc+u°) if Uc > uj o or Uc = Wo; an d faster (~ e -#K+^ ) 2 ) jf trie ener gi es 

are comparable but inequal. We observed similar behavior for the bulk two-point 
amplitude (22) at high energies in (5l. 



3.3 Boundary-boundary amplitude 



Finally, we shall analyze the boundary two-point amplitude. Momentum conservation 
fixes k\ = — k\ = A;" and from the on-shell conditions for low mass excitations — u 2 + 
k 2 = O (a; ) we get uj\ ~ U2 = uj at high energy. The electrostatic two-point partition 
function at equal charges was found only numerically for general r = T2 — T\ in j4]. 
Therefore, we shall calculate the amplitude in the equilibrium configuration. Notice 
that due to symmetry the configuration where the charges lie at antipodal points, 
t~i = T\ + 7T, is always a solution to the equilibrium equations (52): if we set T\ = 
and T2 = 7r the charge distribution is symmetric with respect to the real axis, and 



the imaginary parts in both of the terms of (52) vanish. The total energy is found by 
setting £i = £ 2 in Eq. (18) which yields 



£(2pt,s) = _p( N + 



±F(N) -\F(N + 40 +e log 4£ 



(29) 
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Including the spatial momentum dependence from Eq. (|8j), and by using Eq. (11), we 
find 



A)+2 M 



in (2ttA) 



— 2< % ") exp f-2 £ (2pt ' s) 



N=2iuj,^=—iui 



sinh 27ru; 

After using the on-shell condition — u 2 + k 2 — O (ui°) the amplitude becomes 

smh 2nu 



(30) 



(31) 



The result vanishes for large energies as A ~ e~ . It is in accord with the one sug- 
gested in 21 and also matches with the bulk-boundary amplitude (25 ) asymptotically 
at lo c = oj . 



4 Higher-point amplitudes 

4.1 Bulk-boundary-boundary amplitude 

We can extend our method also to the three point amplitude A 1+2 . As above, we 
place the bulk charge at the origin, where it decouples from the Coulomb gas analysis. 
We then fix the boundary charges at the equilibrium configuration, where they are 



at antipodal points on the circle. The relevant total energy is given by Eq. (18). By 



using Eqs. rt8h and (11), we find 



A / \ W ( 27r ^) ! ^°"'' 1 n2?.fc 11 [ r. cito<£\\ 



(32) 



where the subscripts 1,2 (c) refer to the open strings (closed string), and a = 1,2, c 
in the sums. The 26-momenta can be written as 



k 2 



(wi,A:},0) 
(u 2 ,kl,0) 



(33) 



At high energies mass-shell conditions give u a 



a = 1,2, and therefore 



fcf + k\ ) + ( k x ] > cjj + lo\ + 2wiW 2 cos > {uj\ — Lo 2 y 



(34) 



where <fi is the angle between the vectors k\ and k\. The result for the amplitude may 
be written as 



^1+2(^0^1,^2) 



in (2n\) 



sinh (vr £ a w a ) 

x exp{(wi+w 2 +w c ) 2 [H(ai, a 2 , 0) ± z7r6(o;i, a 2 )] + O (logo;)} 



(35) 
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Here the functions H and are defined as 

1 



ciai, a 2 , 



[F(|l-2ai-2a 2 |) + F(l-2a x ) + F(l-2a 2 )] 



+ - [F(2 ai +2a 2 ) + F(2ai) + F(2a 2 )] 

-2F(a! + a 2 ) + 2F(l-a 1 -a 2 ) - 2 (l-cos0) a!a 2 log2 (36) 



6(«i, a 2 ) 



a i + a l ~ \ (l-2ai-2a 2 ) 2 9 (2a 1 +2a 2 -l) 



(37) 



where F(x) = (x 2 logx)/2 as usual, is the step function, and a a = u a /J2b UJ b- 



Notice that the second inequality in Eq. (34) restricts < a a < 1/2, i.e., neither of 



the open strings can alone carry more than half of the total emitted energy. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 

ari 

Figure 1: The function H for cos0 = 0.6. The curves lie at constant H as indicated by 



the labels, and the striped region is the kinematically allowed one, given by Eq. (38). 



It is crucial that the function H is negative in the physical region for the result to 
make sense: the amplitude must not diverge at high energies. The first inequality in 



Eq. (34) may be written as 



2aio; 2 (1 — cos0) > 2a.\ + 2a 2 — 1 . 



(3* 



We verified numerically that S is negative (or zero) whenever this inequality is sat- 
isfied. Fig. [T] shows the situation for cos0 = 0.6: H is indeed negative in the whole 
kinematically allowed (striped) region. This is a most nontrivial check of our result. 
Notice that H vanishes if cos0 = ±1 and Eq. (38) is saturated. In fact, H vanishes if 



and only if the spatial momenta of all the strings are (asymptotically) parallel, with 
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the understanding that zero is always parallel to any vector. Thus the amplitude 
decays exponentially 



~ e 



for large energies in these configurations, and the decay is even faster 



A 1+2 (u) e ,Wi,W2) ~ exp 



(39) 



(40) 



in other cases. 



4.2 Boundary amplitudes 

We can also give other conjectures on the asymptotics of n-point amplitudes with 
n > 3. We are, however, limited to special kinematic settings, which can be accessed 
by solutions for the symmetric potential problems in |4]. For example, let us consider 
"pairwise" production of four open strings, with k\ 
possibly different energies u a 



ka 



■ k^ , and k 2 



— fcjj, and with 



a = 1,2. We fix the charges symmetrically as 
explained above before Eq. (19). In the case of production to orthogonal directions 
(kf ■ k\ = 0) this configuration also solves the equilibrium equations (52). Proceeding 
as above with £( 4pt ) from Eq. (19), 

in (27rA) 2lK+ " 2) 



^0+4(^1,^2) 



exp [F(2wi) + F(2u 2 ) 



sinh [2ir {u\ + io 2 )\ 
-2F{ut + w 2 ) - 2F(\u 1 - co 2 \) ~ 2 {w\ + w|) log 2] 



(41) 



One can check that the result vanishes asymptotically for any fixed ratio as 
A ~ e -#(^i+^) 2 _ 



5 Conclusions and outlook 

In this paper we used electrostatic techniques to investigate string scattering ampli- 
tudes of D-brane decay in bosonic string theory. In particular, we studied the high 
energy limits of open and closed string emission amplitudes in the half S-brane back- 
ground. We considered pair production of open strings and closed strings. We also 
derived a result for a mixed amplitude with one closed string and a pair of open 
strings, and briefly discussed n-point open string amplitudes with n > 3. 

Overall, our analysis revealed the expected exponential fall-off behavior at high 
energies - the amplitudes decay with sums of the emitted energies in the exponent. 
However, in many cases the decay was found to be even faster, depending on the 
conditions for the spatial momenta. 

An attractive feature of the electrostatic method is that it provides intuitive insight 
into the high-energy behavior of the string amplitudes. It would be worthwhile to 
generalize our investigations to other unstable systems. Some cases to study are: 1) 
full S-brane, which corresponds to a collection of positive and negative unit charges 
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Jl4|, 2) non-BPS Dp-branes for p odd/even in Type IIA/B superstring, corresponding 



to paired Coulomb gases 15 , and 3) inhomogeneous or lightlike decays, possibly 
corresponding to two sets of distinct Coulomb gases. The continuum limit with 
appropriate external bulk or boundary charges in each case would help to find an 
approximate high energy emission amplitude for closed or open strings. 

The high energy closed string pair production amplitude which we obtained is cur- 
rently the only explicit result for this process. There are two remaining puzzles which 
we have so far failed to solve. First, in order to maintain symmetry in exchanging 
the closed strings we had to fix the energies of the closed strings to be equal [5] . This 
requirement is a limitation. It does not arise from the electrostatic approximation - 
in the exact power series expression of the amplitude, each term in the expansion is 
by itself asymmetricH However, it would be natural for the final amplitude to depend 



symmetrically on the energies of the emitted closed strings. 



The second puzzle is associated with open-closed duality 22 . (For discussion 



on the issue, see e.g. 23-26].) A striking mismatch of the duality is the failure to 
connect the closed string IR channel to the open string UV channel. Consider an open 
string annulus diagram, with rolling tachyon background on a decaying p-brane at 
both edges. The two natural ways of cutting the annulus, and the two corresponding 
different kinematical limits, give total amplitudes for open string or closed string 
pair production, with UV region of the open string channel corresponding to the IR 



region in the closed string channel 23 . However, the closed string production rate 
is finite for p < 23, whereas the open string pair production rate in the UV is finite 
for p < 22. In an ongoing work, we have tried to improve the open string analysis by 
including logarithmic corrections to the exponent, using the electrostatic approach, as 



in (31 ). However, the mismatch between the open string and closed string production 
rates seems to become even more pronounced. We have estimated the correction 
numerically and have found the open string pair production rate become UV finite 
for all Dp-branes, with a bound p < 26Q So the question remains, is it valid to think 
of the bulk one-point and the boundary two-point amplitudes as coming from the 
same vacuum open string one-loop amplitude? 

There are some caveats. First, there is no rigorous justification of the analytic 
continuation method from Euclidean to Minkowski signature, proposed in 21 , to 
obtain the exact open string pair production amplitude. Furthermore, it seems to 
be very difficult to extract the IR limit of the open string amplitude in 21 , in 
order to make contact with the closed string UV channel. As far as we know, no 
results in the open string IR channel are known. Second, there are no results for 



amplitudes in superstring theory beyond the bulk-boundary amplitude 18 . Bulk 



one-point amplitudes have been calculated in 27 , and the closed string production 
rate in the IR can be easily extracted to be finite for p < 7. It would be interesting 
to generalize our electrostatics methods to Type II superstring and find out how the 
open string pair production amplitude behaves at high energies. 



6 Notice that each term in the expansion is highly off-shell, so one does not really expect individual 
terms to be manifestly exchange symmetric. 



7 The emitted energy still diverges for p = 25. 
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A Various contractions 

In this appendix we will fill in some gaps between ^ and For ease of reference, 
let us record (|5j, 



N=0 a=l a=n c +l 



AT , / N n c +n a 

at. 



II \ II e^^ J] e^' x ' ( ^M . (42) 



X . . 

2tt , 

i=l \i=l a=l 



To calculate the full contraction, it is useful to include various contributions one 
by one. Denote 



ln c +n 



K nc+no = ( II e^-*>-'*»M . (43) 

Let us first focus just on closed strings, n = 0. It is important to recall that the closed 
strings have mixed boundary conditions, Neumann for parallel ones (/x = 0, . . . ,p) 
and Dirichlet for perpendicular directions (// = p + 1, . . . , 25). This is encoded in the 
Green's functions [61 



(X'»(w a ,w a )X' l '(w b ,w b )) 



Neumann 
Dirichlet . 

(44) 

These yield (the singular self-contractions are dropped) 



r -i^(i g 


\w a - w b \ 


2 + log 


11 


- w a w b \ 


\ -|^(log 


\w a - w b \ 


2 -log 


ll 


- w a w b \ 



K nc+ Q= |[ \1 - W a W b \ ka ' k b k "' k b\w a - Wb \ ka ' k b +k *' k b Y[ |1 - W a W a \^ ka)2 2^) 2 . 

l<a<b<n c a=l 

(45) 

Now we wish to take into account open strings, i.e., n Q ^ 0. They only couple 
to the parallel parts of the fields and have the Neumann boundary conditions. The 
contribution is thus easily accounted for: 



n c n 



K nc+no = K nc+0 n K-^n' fc TI 11 K-^n^. (46) 



n c +l<a<b<n a=l b =n c +l 
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Notice that there is a factor of 2 relative to bulk-bulk case in the exponents, since 



the two terms in (44) add up. 



Finally, we wish to include the contribution from the boundary deformation. Since 
the deformation only involves the field X°, we get the contribution only from the 
temporal direction: 



N n c +n 
no 



n 



i i ^ 



a=l 



l<i<j<N 



N n c +n 

n n k 

i=l a=l 



Jvt% I 2^ a 



(47) 



Inserting (47) to (42), with the expression for K nc+rio (46) and for K nc+ o in ( |45 
yields Q. 



B Equilibrium conditions 

We shall look for the (global) equilibrium configuration, which is the electrostatic 
configuration in the Coulomb gas picture. Let us start with the boundary n-point 
amplitude A 0+n and set all spatial momenta are zero, k a = 0. Since the external 
charges lie on the unit circle, w a = w~ . Therefore, we may extend the Coulomb gas 
partition function to an analytic function of w a : 



Z{N) 



1 



N 



n£ n 

2=1 l<i<j<N 



I g^i I ^ 



w a ) + \og(e 



x ex p { Yl ^ [ log ^ ei 

i,a 

+ [log{w a - w b ) + \ogiw~ 1 - w^ 1 )] I 



(4J 



We find that the saddle point equations can be written as 



O=i0„,logZ(JV)=& 



E 



+ N 



W n 



W r 



W f 



2w n 



(49) 



where the expectation value is defined by Z(N). Notice that (49) is basically (the 
expectation value of the conjugate of) the electric force felt by the particle at w — w c . 
The first term is due to the self-interactions of the external charges at w = w a , and 
the second term is the expectation value of force due to the unit charges created by 
the tachyon profile. This suggests that the solution of the equations for all w c is the 
equilibrium configuration. 



The last term in (49) is an electric force due to a special charge at w = 0. The 



origin of this term is understood as follows. Note that (48) is real (for real £ a ) when 
all w a are on the unit circle. Hence the complex derivative with respect to any w c 
must return an tangential force, i.e., oc iw~ l where the proportionality constant is 
real. The radial force equation is automatically satisfied when all \w a \ = 1. In Eq. 
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(49), this is explicitly realized by the additional charge at the origin, which cancels 



the radial pressure due to the interactions of the charged particles. 



We may verify these arguments explicitly by splitting (49) into radial and tangen- 
tial components. We write 



W 



w. 



2 d Wc \ogZ(N) + ^d Wc \og Z(N) 



N 



E 



a=£c 



W n - W, 



+ 



IWc 

2 

i 

N 



d Wc \og Z(N) 



W a - W, 

iw c 



+ 



e ttl _ w 



+ 



w c 



+ 



N 



d w log Z{N) 



W n 



Wr 



W„ 



w c 



+ i 



Wr 



W r 



W r 



Wr 



(50) 



where the former expression is the radial force and the latter one is tangential. Since 
both w a and w c lie on the unit circle, 



W r 



+ 



1 



W n 



\W r 



1 • 



(51) 



w a — w c w a - w c 

i.e., the radial electric field at w c due to a particle at w a is independent of both 



w a and w c . Hence one sees immediately that the radial component in (50) vanishes 



identically, so the equilibrium configuration is fixed by the tangential equation. If one 
uses rotational symmetry to fix w c = 1 the tangential equation becomes 



1 T 1 
-)^£Jm - + 

N w a -l 



Im 



. 



(52) 



Above derivation was done for fixed N. When the partition function is summed 
over N we expect that the final equilibrium equations are found by continuing an- 
alytically to iV = — ^2 a £, a , m analogue with the partition function. In the text we 
shall apply the equations only to such cases where the saddle point configuration is 
independent of N. 

The case of nonzero k a is also interesting. As is easy to see from (JsJ) , for the 
boundary amplitude this means replacing £ a by an effective charge £a = ^ a + k\-kl/^ c 
in the above formulas. Extension to bulk charges is simple as well. In this case w a and 
w a can be taken to be independent, and the electric force obtained by differentiation 
has two nontrivial components. 



References 

[1] D. J. Gross and P. F. Mende, Phys. Lett. B 197 (1987) 129. 
[2] D. J. Gross and P. F. Mende, Nucl. Phys. B 303, 407 (1988). 
[3] D. J. Gross and J. L. Manes, Nucl. Phys. B 326, 73 (1989). 



16 



[4] N. Jokela, M. Jarvinen and E. Keski-Vakkuri, [arXiv:1003.3663 [hep-th]. 



[5] N. Joke la, M. Ja rvinen and E. Keski-Vakkuri, Phys. Rev. D 80 (2009) 126010 
[arXiv:0911.0339| [hep-th]]. 

[6] V. Balasubramanian, E. Keski-Vakkuri, P. Kraus and A. Naqvi, Commun. Math. 
Phys. 257 (2005) 363 |arXiv:hep-th/0404039| . 



[7] S . R. Das and S. D. Mathur, Nucl. Phys. B 482, 153 (1996) |arXiv:hep- 
|th/9607149| . 

N. Jokela, M. Ja rvinen and E. Keski-Vakkuri, Phys. Rev. D 79 (2009) 086013 
[arXiv:0806.149T| [hep-th]]. 

[9] A. Sen, JHEP 0204 (2002) 048 [arXiv:hep-th/0203211] . 



[10] F. Larsen, A. Naqvi and S. Terashima, JHEP 0302 (2003) 039 |arXiv:hep- 
|th/0212248|. 



[11] N. D. Lambert, H. Liu and J. M. Maldacena, JHEP 0703, 014 (2007) |arXiv:hep- 
|th/0303139| . 

[12] V. Balasubramanian, N. Jokela, E. Keski-Vakkuri and J. Majumder, Phys. Rev. 
D 75 (2007) 063515 [arXiv:hep-th/0612090] . 

[13] N. Jokela, M. Jarvinen, E. Keski-Vakkuri and J. Majumder, J. Phys. A 41 (2008) 



015402 |arXiv:0705.1916| [hep-th]]. 
[14] N. Jokela, E. Keski-Vakkuri and J. Majumder, Phys. Rev. D 77 (2008) 023523 



|arXiv:0709.1318 [hep-th]] 



[15] J. A. Hutasoit and N. Jokela, Phys. Rev. D 77 (2008) 023521 |arXiv:0709.1319 
[hep-th]]. 

[16] N. Jokela, M. Ja rvinen and E. Keski-Vakkuri, Phys. Rev. D 79 (2009) 106005 
|arXiv:0901.3368| [hep-th]]. 

[17] A. B. Zamolodchikov a nd A. B. Zamolodchikov, Nucl. Phys. B 477, 577 (1996) 
[arXiv:hep-th/9506136] . 

[18] N. Jokela, E. Keski-Vakkuri and J. Majumder, Phys. Rev. D 73 (2006) 046007 
[arXiv:hep-th/0510205] . 

[19] N. Jokela, M. Jarvinen and E. Keski-Vakkuri, J. Phys. A 41 (2008) 145003 



|arXiv:0712. 4371 [cond-mat.stat-mech]]. 



[20] M. E. Fisher and R. E. Hartwig, Adv. Chem. Phys. 15 (1968) 333-353; H. Widom, 
Amer. J. Math. 95 (1973) 333-383; E. Basor, Transactions of the American 
Mathematical Society 239 (1978) 33-65. 



17 



[21] M. Gutperle and A. Strominger, Phys. Rev. D 67 (2003) 126002 |arXiv:hep- 
|th/0301038| . 

[22] A. Sen, Phys. Rev. Lett. 91 (2003) 181601 | arXiv:hep-th/0306137| . 

[23] J. L. Karczmarek, H. Liu, J. M. Maldacena and A. Strominger, JHEP 0311 
(2003) 042 |arXiv:hep-th/0306132| . 

[24] A. Sen, Int. J. Mod. Phys. A 20 (2005) 5513 | |arXiv:hep-th/0410103] . 



[25] Y. Nakayama, S. J. Rey and Y. Sugawara, JHEP 0608 (2006) 014 |arXiv:hep- 
|th/0605013| . 

[26] Y. Song, JHEP 1007 (2010) 030 [arXiv:1003.0230| [hep-th]]. 
[27] J. Shelton, JHEP 0501 (2005) 037 |arXiv:hep-th/0411040| . 



18 



